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Abstract 



^1 We compare two methods of analyzing the finite-temperature electroweak phase transi- 

tion in the minimal supersymmetric standard model: the traditional effective potential (EP) 
■ approach, and the more recently advocated procedure of dimensional reduction (DR). The 
' latter tries to avoid the infrared instabilities of the former by matching the full theory to 
an effective theory that has been studied on the lattice. We point out a limitation of DR 
that caused a large apparent disagreement with the effective potential results in our previous 
work. We also incorporate wave function renormalization into the EP, which is shown to 
decrease the strength of the phase transition. In the regions of parameter space where both 
methods are expected to be valid, they give similar results, except that the EP is significantly 
more restrictive than DR for the range of baryogenesis-allowed values of tan/5, mh, the crit- 
ical temperature, and the up-squark mass parameter mjj. In contrast, the DR results are 
consistent with 2 <, tan /? ^ 4, nih < 80 GeV, and mu sufficiently large to have universality 
of the squark soft-breaking masses at the GUT scale, in a small region of parameter space. 
We suggest that the differences between DR and EP are due to higher-order perturbative 
corrections rather than infrared effects. 



1 Introduction 



Recently a new method has been proposed and exploited to try to improve the reliability 
of analyzing first-order phase transitions in finite-temperature gauge theories, called dimen- 
sional reduction (DR) In the traditional effective potential (EP) approach, a problem 
is posed by the transverse gauge bosons which are very light, and thus can lead to infrared 
instabilities in the effective theory. Their effects should in principle be treated nonpertur- 
batively to reliably compute the strength of the transition. Dimensional reduction works 
around this by integrating out all the heavy modes in a given theory to obtain an effective 
theory of the problematic light modes, which is then studied on the lattice . This program 
was carried out in refs. 0] for the standard model. However, the nonperturbative lattice 
results of ref. are applicable to other theories as well, since the only requirement is that 
the effective theory after integrating out the heavy fields has the same degrees of freedom as 
does the standard model, namely the light Higgs doublet and the transverse gauge bosons. 

One of the most interesting applications of DR is to supersymmetric (SUSY) models, 
since these have the possibility of producing the baryon asymmetry of the universe at the 
electroweak phase transition 0. In order to preserve the baryons thus created, the phase 
transition must be strongly enough first order, meaning in this context that the VEV of the 
Higgs fields must be sufficiently large at the critical temperature: 



Otherwise anomalous baryon- violating interactions mediated by electroweak sphalerons will 
be too fast in the broken phase after the transition, and quickly wash out the baryon asym- 
metry that was created. The range of values on the right hand side of (|l|) comes from an 
estimate of the uncertainty in the sphaleron rate at two loops, and in dynamical details like 
the amount of reheating in the phase transition 0]. In the DR approach, condition (|l|) is 
replaced by the requirement that the ratio of the Higgs quartic coupling to the squared gauge 
coupling must be small enough in the effective three-dimensional theory. 



The two quantities are related to each other by lattice measurements of v(Tc)/Tc in ref. [0; 
we have used their results to obtain the fit 



v{T,)/T,^{vl + viy/yT,> 1-1.5 



(1) 



Xc = h/gl < 0.026 - 0.044. 



(2) 



v{Tc)/Tc ~ 1.087^2 - 2.916w; + 2.911; w = xjOM, 



(3) 



which is valid in the interval corresponding to (P and (|^). 
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References have studied the phase transition in the minimal supersymmetric stan- 

dard model (MSSM) using DR, identifying the regions in the SUSY parameter space con- 
sistent with condition (Q). One of the purposes of this letter is to explain the differences 
between our results |Q and the others In the process we shall elucidate a possible 

shortcoming of DR: even moderately heavy squarks can be too light to reliably integrate out 
if the superrenormalizable cubic scalar couplings of the MSSM become too large. 

More generally we are interested in determining the extent to which the results of dimen- 



sional reduction differ from those of the effective potential (EP) approach P]-||Tj]- We have 
therefore undertaken a comparison of the two methods in the MSSM, exploring as broad a 
region of parameter space as possible. We find that they give results which are in reasonable 
qualitative agreement. A further goal is to characterize which regions of MSSM parameter 
space are suitable for electroweak baryogenesis. We have thus recomputed the distributions 
of MSSM parameters, as well as some observables, that are baryogenesis-compatible. Where 
comparison is possible, our results agree well with some of the previous DR studies that 
have focused on more specific regions of parameter space, and less well with others 0. 

The dimensional reduction procedure employed here is essentially the same as we used 
in ref. 0] (hereafter called CK), which includes the one- loop corrections proportional to the 
top and bottom yukawa couplings, ytfi. One must integrate out the third family quarks and 
squarks at zero and at finite temperature to find the effective finite-T lagrangian and to 
relate its couplings to physical quantities. The combination Xs/g^ in eq. (Q) can thus be 
expressed a function of squark masses and mixing angles, Higgs boson masses and tan/5, 
allowing one to infer from the inequality (|^) which are the baryogenesis-compatible regions 
in the space of physical parameters. 



2 Dimensional reduction: subleading corrections 

Because Xc needs to be small for purposes of baryogenesis, we are interested in cases where 
there are large cancellations between the tree- level and one-loop contributions to A3. For 
quantitative accuracy it therefore appears important to include some formally subleading 
contributions in the finite-T effective lagrangian, which were omitted in CK. Among the 
potentially most important such terms are those due to thermal loops of gauge and Higgs 
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bosons, which shift the Higgs boson mass terms by an amount 

= + ^Tm^ + \H,n (4) 



= cos22a - TTT^^::— , (6) 



where 

These additional terms lower the estimate of the critical temperature from that found in 
CK; they alone can increase the final value of the parameter Xc by ~ 0.008. However, there 
is also another, direct contribution to Xc from the same particles, 

gHn2 2 3g'T 

— cos 2a . 

327r2 USirmAo 

where the angle a defines the direction of the phase transition in the {Hi, H2)-pla.ne. Both 
in (^) and (^) the first term comes from superheavy (nonzero Matsubara frequency) and 
the second from heavy scale (zero Matsubara frequency) integration. The correction (|^) is 
~ —0.003, partly cancelling the effect on Xc from the thermal self-energies in (^; we argue 
that this combined effect of ~ 0.005 is the typical scale of dominant (^^-corrections. 
Taking into account these improvements, our full result for be written as 

Xr = — — COS 2a 




+ ( — sin a ^-^ — cos 2q; sin a 



( C(3) 



T 



/i^ cos^a — Af sin^a) cos 2a 



+ (t ^ b; cos a ^ sin a) , 

+ 6x^^ + 6x" + 6xg (7) 

Here M^t is the sum of the Debye masses left- and right-handed top squarks [^, 

Mn,=mi^+mi^, (8) 

and the thermal corrections are computed assuming the gauginos and higgsinos are decou- 
pled. The first three lines in (0) are, respectively, the contributions from tree-level, leading 
superheavy scale, and the superheavy and heavy scale wave function renormalizations. The 
fourth line shows how to obtain the effect of the bottom sector from that of the top (no- 
tice that the interchange of sin a and cos a implies cos 2a — cos 2a). The last line gives 
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the next-to- leading superheavy scale shift Sx^^ , the much larger term Sx^ from the same 
Feynman diagrams in the heavy scale integration, and the gauge boson term from (^. The 
additional function 6x^^ is given by 

1/3 

9 ^ ^ a=Q,U 

3 3 

- i^yt{Slffi^{^Q^^Q^^u,mu) - -g'^ytF^^{mQ,mu) cos 2a 



+ {t ^ b] U ^ D; sina ^ cos a) 



fl/f (^Q, rnu) - ^F!"{mu, m^)) 



+ ^(mQ, ttid) + 2Ff ^(m^, mg)) j cos 2a (9) 

with F^^{ma,mh) = ?,ir)?af^^ {nia.nia) - S^^f^^ {ma,ma,mb), S*^ = (yucosa + A sina)^ 
The high-temperature expansions of the n-point loop integrals f^^ (following the notation 
of [|) are 



1287r4T2 1287r6 
C(5) 

/f^(m«,ma,mb,mb) ~ ^^^^^6^4 (10) 

to leading order in m? /T"^. The contribution from the heavy scale, 5x^ , has the same form 
as 6x^^; one merely replaces the masses mq uD by the corresponding Debye masses (see ref. 
@]) and the integrals (|l3) by their heavy scale counterparts, which are exactly given by 

T 

f 2 {ma, ma) = 



■inrua 



h {ma,ma,mb) ,^^^3 

T 



Snniaima + m;,)^ 

f^{ma,ma,mh,mb) = f ■ r^r. (11) 



^nniamhima + mb)^' 

In CK we omitted 6x^^ and the corresponding corrections to the wave function renormal- 
ization, since they are formally subleading. The combined effect of all the new terms we 
consider here is an increase in Xc of varying size, but never larger than 0.015. 
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Figure 1: A contribution to the effective quartic coupling of Higgs bosons involving the 
interactions of eq. (|13), and generating the term ( p!3[ ) in the effective Lagrangian. At finite 
temperature it is infrared divergent if the masses of the squarks in the loop vanish. 



3 Constraint on superrenormalizable couplings 

Despite the significance of the subleading contributions which we now take into account, the 
most important difference between this work and CK is that we have imposed a consistency 
condition on the size of the fi and f, parameters of the MSSM, the necessity of which 
was not recognized previously. These couplings are distinguished by the fact that they 
appear in superrenormalizable interactions among the squark and Higgs bosons in the MSSM 
Lagrangian: 

>Cmssm = yt ilijjHi + AtH2)iR + bl{fiH2 + AbHi)bR + ■■■ (12) 

Because of their dimensionality, one must be especially vigilant against the breakdown of 
the heavy scale perturbative expansion due to large infrared contributions to the three- 
dimensional finite-T effective Lagrangian (where the effective quartic coupling has dimensions 
of T), of order 

SX3-T^^\f^H, + AH2\\ (13) 

coming from diagrams like figure 1. In the present study we have imposed an upper limit 
on the ratio of fi and At^i, to M^^ ^ to safeguard against such a breakdown. Precisely, we 
required that 



b\ 



< 1; and \fi\, < 500 GeV. (14) 



t,b 



The precise values chosen for the bounds (|Hp are somewhat arbitrary. However the region 
of MSSM parameter space that satisfies the sphaleron bound but violates ( P^^D appears to 
be an isolated "island" which, once excluded by ([14D , is not further reduced by making ( |14|) 
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more stringent. These new consistency conditions make a dramatic change in the regions 
of MSSM parameter space which are most hkely to give a strong phase transition, and are 
the main reason for the differences between the present results and those in CK. There 
our baryogenesis-allowed points were dominated by surprisingly high critical temperatures 
near 300 GeV. We have subsequently investigated such large-Tc points using the EP and 
found that for them, v/T falls below 1 very rapidly with a slight increase in T. Since the 
present estimate of (the temperature where the second derivative of the potential at zero 
field has vanishing determinant) is already an underestimate of the true (where there are 
degenerate minima of the potential), such behavior means that the transition is quite unlikely 
to be strongly enough first order for baryogenesis. We find a clear separation between these 
spurious solutions to the sphaleron bound, and the acceptable ones with near 100 GeV, 
on which we focus henceforth. 

The conditions ([T^) also constitute an important difference between DR and EP. In EP 



the constraints (14) are not needed, unless one made an expansion of the field-dependent 
mass eigenvalues of the squarks in powers of At and /i. But since EP uses the exact field- 
dependent tree-level squark masses, it resums what would be an infinite series of nonrenor- 
malizable contributions to the effective theory.[] In DR on the other hand, one must match 
the parameters of the effective theory onto the renormalizable Lagrangian that has been 
studied on the lattice. We therefore have no choice but to truncate the effective theory at 
the order of the renormalizable operators. To improve upon this, it will be necessary to 
simulate the full MSSM on the lattice, but in the meantime EP is the only method which 
has a hope of reliably determining the strength of the phase transition for values of /i and 
At where (|lj) does not hold. 



4 The effective potential 

To compare DR to the effective potential approach, we constructed the EP at the one-loop, 

ring-improved level, including contributions of the virtual standard model particles plus 

the top and bottom quarks and squarks, and the Higgs bosons.0 This extends the work 

described in ref. where the bottom sector and Higgs bosons were omitted. We find that 

^See the caveat discussed below about resumming the corresponding thermal corrections to the squark 
mass matrices 

^The Higgs bosons require special treatment in the EP because some of them can have negative for 
small values of Hi. We dealt with this by setting the contributions to the EP to zero whenever m? < 0. 
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the inclusion of the bottom squarks shifts the critical temperature significantly (10 — 20%), 
but has a small effect on f/T, for experimentally allowed values of the lightest Higgs boson 
mass.Q After successfully reproducing the results of ref. 0, we further improved the EP by 
supplementing it with the same wave function renormalization as we computed in the DR 
approach. Usually wave function renormalization is ignored in the EP, and the one-loop, 
ring-improved effective lagrangian for the Higgs fields is taken to be 

£efT = El^.^^l'-^tree(^0 + ^StrM^(if,)(^ln^^^^-^^ 

+ StrT|^ln (lie-V^-Wf/,)^ _ _Z_Tr(M|,(/f„T)-M3(ff,)). (15) 

Here Str denotes the supertrace, ± is +(— ) for fermions (bosons), Tr is the trace over 
bosons only, itld is the thermally corrected Debye mass, and Q is the renormalization scale. 
However it is more accurate to also include the renormalization of the kinetic term so that 
Y,i l-D^-f^jP becomes Y^ij ZijiD^HiY D^Hj. After rescaling the fields to canonical form and 
ignoring effects of two-loop order, this amounts to making the replacement 

j 

One reason for making this improvement is to try to minimize any possible sources of discrep- 
ancies between DR and EP. We find that including wave function renormalization reduces 
the strength of the phase transition noticeably, though not drastically. The comparison is 
shown in figure 2, where we have plotted contours of v{Tc)/Tc in the plane of tan (3 and 
niA (the pseudoscalar Higgs boson mass), for nit = 170 GeV, ttiq = 280 GeV, mu = 0, 
At = iJ = 0. 

To include the effects of the Higgs bosons in the EP, we used the following for the 
field- dependent parts of the mass matrices for the CP-even, CP-odd and charged bosons, 
respectively: 



M'hW = Mi(0) + ^^ 2^' + '^' ° 





4 I 3y^ + g' 




or 



Hf + Hi AH1H2 
AH,H2 Hi + Hi 




■^We explored very large values of tan /? ~ 45 where the bottom quark Yukawa coupling would be relevant, 
but found no cases with v/Tc > 1 except those suffering from the same problem as the large Tc points we 
rejected in connection with (|l4|). 
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without w.f.r. with w.f.r. 




100 200 300 400 500 '"^ 100 200 300 400 500 '"^ 

Figure 2: Comparison of the usual effective potential results for the phase transition in the 
MSSM, to the EP supplemented by wave function renormalization. Heavy convex lines: 
contours of constant v/T at the critical temperature, in the tan j3-m a plane. Heavy concave 
line: contour for ruh = 65 GeV. Light concave lines: contours for other values of rrih- Fig. 
2a is without wave function renormalization, and 2b is with w.f.r. The values of the other 
MSSM parameters are nit = 170 GeV, mq = 280 GeV, mu = 0, At = = 0. 



5 The high-temperature expansion 

To obtain the above results for the EP, we used the exact expressions for the T-dependent 
part, without expanding in masses over temperature. In contrast, DR is explicitly a high 
temperature limit, requiring us to impose the important restriction 

m^^/T < TT (18) 

on the parameters considered, where rriq^ is the thermal (Debye) mass of any of the squarks 
included in the loop corrections, qx = iL,R,bL,R 0- One might therefore naively conclude 
that EP is valid for a larger range of parameters than is DR. However we wish to argue 
that one should also make the same restriction (0) to obtain reliable results from the EP, 
because the ring resummation assumes that all the modes in the "soft" loop (heavy scale) can 
be taken approximately massless, compared to the modes in the "hard" loops (superheavy 
scale). Thus the ring-corrected EP is strictly speaking, and somewhat contrary to common 
wisdom, valid only in high-T limit. Since the EP has the qualitatively correct decoupling 
limit, when one employs the nonexpanded integral expressions for loop-corrections however, 
one might hope that the infinite resummation retains some of the qualitative physical features 
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of the exact theory, even outside its region of strict vahdity. Nevertheless, in what follows 

we will impose (|T8|) equally on DR and EP. 

At the same time however, we cannot allow any of the particles we are integrating out to 

become arbitrarily light compared to T, because this would give rise to infrared divergences, 

as discussed above, and necessitate a new lattice study with the light fields included in the 

low-energy lagrangian. Thus we imposed the further constraint 

T ~ ~ 
< 1; qx = tL,R, hL,R. (19) 

This ratio appears as an explicit expansion parameter when integrating out the heavy-scale 
degrees of freedom; it is the exact analog of the expansion parameter gT /Mw in the gauge 
sector, which at small M\y ~ gcf) induced us to use DR rather than EP in the first place. We 
found that EP is quite sensitive to the exact value taken as the upper limit in (|r^): using 
a value of 1 removes 80% of otherwise acceptable parameter sets, whereas using 1.2 leaves 
essentially all of them. In the present work we imposed (19) only on DR, not on EP. 



6 Monte Carlo search of the parameter space 

We undertook a Monte Carlo sweep of the MSSM using both methods, DR and EP, to find 
those parameters allowed by the baryogenesis requirement (|l|) or (|]) 

< 0.044 ^ v{T^)/T^ > 1.0, (20) 

and also the bound on the p parameter, which we took to be Ap < 0.011 as the contribution 
from the third generation quarks and squarks. It can be seen below that the scarcity of 
solutions to ( ^OD depends strongly on small changes (~ 0.001) in the value of this upper 
bound, so that future improvements on the p parameter constraint might severely limit 
the possibilities for electroweak baryogenesis in the MSSM. As in CK, the randomly-varied 
independent parameters were tan/3, rrit, the pseudoscalar mass uia, the p parameter, and 
the soft SUSY-breaking parameters mg, mu, mci, At and Ab. The derived quantities are 
the physical masses of the squarks and lightest Higgs boson, the critical temperature Tc, 
Ap, and v/T or Xc- To compare results, we have converted Xc into the equivalent v/T value 
using eq. (^. To determine v/T with the EP, we numerically solved for the global minimum 
of the potential, but found that for the parameters of interest this was never more than 
1% from the value obtained by minimizing the one-dimensional slice of the potential in the 
symmetry-breaking direction determined at the origin. 
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The resulting distributions of parameters are shown in figure 3, where one sees reason- 
able qualitative agreement between the two methods. However there are several noticeable 
differences. The maximum allowed up-squark mass parameter^ is 160 GeV in DR and only 
90 GeV in the EP. The largest allowed values of m^o, the mass of the lightest Higgs boson, 
are 84 GeV and 70 GeV respectively. In EP we find a somewhat lower critical temperature 
(99 GeV versus 92 GeV), and the distribution of Tc values is also much narrower there than 
in DR. The narrowness of the tan/? distribution in EP relative to DR has the anticipated 
correlation with that of rriho since the latter increases in the MSSM with larger tan (3. Finally 
one sees also that DR has a broader distribution of v/T than does EP. The fact that DR 
has less difficulty than EP to satisfy the sphaleron constraint is evident in all the parameter 
distributions where a difference is discernible. 

As for why DR and EP do not agree exactly, we note that the perturbative expansions 

are somewhat different in the two approaches. In particular the ring- improvement in EP is 

effectively a nonanalytic resummation with regard to the scalar field dependent terms. The 

DR counterpart of this is the heavy scale integration with the SH-scale corrected particle 

(Debye) masses, which is only done to the order of renormalizable terms in the lagrangian. 

Therefore the two methods differ by contributions that are of two-loop order, and also by 

nonrenormalizable terms; we believe these are the source of the discrepancies. It is already 

known that in EP the two-loop contributions may considerably strengthen the order of the 

transition T^\, extending the allowed values in tan/3 up to ~ 4, just as we are now 

finding in DR but at only one loop. There may be significant higher-order corrections to 

DR as well. Although we imposed the constraint (|T9t), to ensure the overall convergence of 

the heavy-scale perturbtion expansion, the two-loop correction may be sizeable when mu is 

small. Indeed one may roughly estimate that 5xc,2-ioop ~ (z/t /47r)^(T/m(-^)^ <, 0.01. 

^ We excluded negative values of mfj in order to avoid color-breaking minima, although it is possible to 
be less restrictive However we did include the empirical constraint + < 7.5(m? + m~ ) from 



reference iTS , since color-breaking minima can occur even when mfj is positive 



+ 
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1.65 1.75 1.85 1 2 3 4 5 10 1 2 3 4 1 

tan^ Mq 




Figure 3: Distributions of MSSM parameters satisfying the baryogenesis constraint in the 
dimensional reduction approach (DR) and in the effective potential approach (EP). Masses 
are in units of 100 GeV. 



7 How small must mu be? 

One might wonder how much fine tuning of parameters is needed to make the phase transition 
strongly first order. One simple assumption is that all the soft-breaking masses are equal 
at the GUT scale, rriQ = mu = mo, and therefore only differ from each other at the weak 
scale by logarithmic corrections. Using the renormalization group equations to run the 



soft-breaking masses down from their universal GUT-scale value |T6[, and keeping only the 



corrections due to the large top quark Yukawa coupling, one can estimate that at the weak 
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m / m 
U Q 



m / m 
u Q 




0.45 0.9 0.3 0.6 

Figure 4: Distributions for mu/mq in DR and EP. 



scale 

mu I Mgut 



0.7-0.8, (21) 

rtlQ V Mweak / 

using yl = 1.25 (corresponding to tan/3 = 2) and Mgut = 10"'^'^ — 10"'^^ GeV. 

It is encouraging that this goes in the right direction to be consistent with electroweak 
baryogenesis, since the first order phase transition requires small values of mjj. However the 
values of mu/rriQ derived from the EP tend to be smaller than this: the maximum is 0.8, and 
values greater than 0.4 are unlikely. In DR, mfz/mg has a somewhat less restricted range, 
going up to 1.6. Although values as large as in (0) are relatively unlikely, at least from the 
point of view of the parameter space for the MSSM unconstrained by GUT relations, they 
are nevertheless possible. The distributions for mu/mQ in DR and EP are shown in figure 4 
(the full extent of the tails of these distributions is not shown). Out of 12,500 accepted points 
in the DR Monte Carlo, about 200 satisfy (0). Among these points, the two top squarks 
are typically split by less than 80 GeV with 190 < m^^ < 280 GeV , the pseudoscalar Higgs 
boson is heavier than m^o > 180 GeV, and the right-handed top squark mass parameter lies 
in the range 60 GeV < mu < 150 GeV. For these points fi and At can be several hundred 
GeV but they conspire to give a small value of fi + At tan f3 so that the top squark splitting 
is small. This region can have a Higgs boson as heavy as rrih = 79 GeV. 



8 Correlations between parameters 

In our previous work (CK), there was a strong correlation between the allowed values of 
tan/3 and tua, including points with arbitrarily large values of tan/3. These points were 
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Figure 5: (a) Baryogenesis-allowed region of MSSM parameter space in DR, projected on 
(a) the plane of m^o versus tan (3 and (b) in mu versus tan j3. 



all associated with large T^, or equivalently with all SUSY breaking parameters at the TeV 
scale, as discussed above. Constraints (|T1| and O) remove all these points from our sample 



and we find the weaker correlation shown in fig. 5a. The correlation between large tan j3 and 
small mu in figure 5b suggests that the failure of EP to allow large tan (3 may be due to the 
effects of having a very small value of uiu. 

Previous studies of the phase transition in the MSSM have emphasized that large values 
of the IX and At parameters have the effect of weakening the transition. While this may be 
true while holding all other parameters fixed, if they are instead allowed to vary, one can 
still find values where the transition is strong even for large /z and At. Despite the restriction 



(0) on /i and At^b that was used to obtain the present results, fig. 6a shows that a large 
fraction of the fi ~ At plane is nevertheless represented. This is worth emphasizing because 
in the MSSM the CP-violating phases which are needed for electroweak baryogenesis are in 
precisely these parameters. Therefore the regions with larger rather than smaller values of 
or \At\ are the more interesting ones. 



It was mentioned above that the p parameter places a stringent constraint on baryoge- 
nesis in the MSSM. One finds that it is strongly correlated with various other parameters; 
in figure 6b we show how in DR it appears to be a limiting factor on how large the light- 
est Higgs boson mass rrih can be. The shape of the correlation between Ap and tan/3 
looks quite similar, as one might expect from the tree- level formula m\ = ^{m\ + m\ — 
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Figure 6: Baryogenesis-allowed region of MSSM parameter space in DR, projected on the 
plane of (a) /i versus At, and subject to the conditions ( p!^p!9D , and (b) on the plane of Ap 
versus nih. 



^J {m\ + m|)^ — Am\m\ cos^2/5): it vanishes at tan /3 = 1 and increases for larger tan/3. 

9 Comparison with other DR results 

In addition to comparing our DR results with the EP, we have also compared to other 
pubhshed DR results. We have analytically checked that our results agree with those of 
ref. 1^, up to different ways of handling the renormalization procedures (which are two- loop 
effects), our keeping the finite g' effects in the loops, and some numerically small additional 
terms included in ref. Numerically our results agree typically to within less than 0.01 in 

More recently, ref. gave results which indicated larger values of the critical ratio Xc 
than did ours, hence a weaker phase transition, for given MSSM parameters. The authors of 
ref. identified the following possible source for the discrepancy: the direction of symmetry 
breaking in the plane of the two Higgs fields is very sensitive to the critical temperature Tc for 
small pseudoscalar mass m^; a large change in this mixing angle can on the other hand cause 
large variations in Xc- We have confirmed this expectation by observing that the tree-level 
value of tan/3 (used by ref. 0) and the one- loop value (used by us) differ significantly for 
small values of mA- Indeed, we find the region of tua < 100 GeV to be that where the results 
of all three groups agree least well. Fortunately this region is not particularly essential for 
baryogenesis, nor for our present results. 
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Full result 
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CK+leading g corr. 

(Full-g') corrections 
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Figure 7: Critical ratio function of niA for tan/? = 1.75, showing the dependence 

on corrections of order g^. Other parameters are mq = 265 Gev and mu = 100 GeV, 



The above observation does not fully explain our disagreement with ref. [0 in the large 
rriA region. In fig. 7 we show Xc as a function of rriA at tan/3 = 1.75, for comparison with 
fig. 2 of For large our asymptotic value of Xc = 0.03 is smaller than theirs by 0.02. 
Of this discrepancy, 0.004 is due to the different definitions of tan/5. We have also studied 
the effects of loop contributions of order to the Higgs boson mass parameters, neglected 
by us but included by the authors of [0], and find that they can give a further increase of 
0.005. We are unable to account for the remaining difference of 0.01. 

One should bear in mind that the natural scale of the leading terms contributing to Xc 
are of the order of ~ 0.2 and the accepted solutions always correspond to a large degree of 
cancellation between these leading contributions. From this perspective we thus have in fact 
a rather good numerical agreement with a relative accuracy of about 10 per cent with the 
other published works in the field. Fig. 7 gives an indication of the sensitivity of Xc to the 
various approximations: omitting all g"^ corrections, as in CK, putting in the leading g^ terms 
but not g'"^, and including all the corrections described in this work. The difference between 
our present work and CK comes dominantly from the neglected squark sector correction 
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Sx^^, equation (H). 

There are several effects that make it difficult to achieve a better agreement at this 
level of approximation. For example we find that differences in the definition of heavy scale 
squark Debye masses introduce changes of order 0.005 in Xc- For both DR and EP we have 
consistently used the leading order approximations for Debye masses, whereas one could 
include also the next-to-leading corrections. In the absence of the T = renormalization of 
the squark sector this would introduce dependence on the renormalization scale p , so to get 
more accurate results the complete renormalization of squark sector at one loop would be 
required. Also, while the other authors were working in the approximation g' = for the loop 
corrections, we did not; this too makes differences of order ~ 0.005 in Xc- Finally, different 
renormalization procedures cause the definitions of ruh to differ by terms of two-loop order. 
If one uses the physical Higgs boson mass as an input instead of tan (3, getting an accuracy 
of 0.005 in Xc would require using the two-loop computation for rrih, since we have checked 
that a difference of 2-3 GeV in rrih changes the corresponding value of Xc by 0.005. 

10 Conclusions 

In summary, we have compared the one- loop dimensional reduction and the effective po- 
tential approaches throughout the parameter space of the MSSM, and found that they give 
qualitatively similar results for the strength of the electroweak phase transition, although 
they differ in certain quantitative respects, which should be further explored. We showed 
that wave function renormalization has a noticeable weakening effect on the phase transition 
when incorporated into the EP. One potentially interesting difference is that DR appears 
to allow a smaller hierarchy between the soft-supersymmetry-breaking top squark mass pa- 
rameters mu and ttlq, which could make SUSY electroweak baryogenesis compatible with 
universality at the GUT scale, mu = mg. The top squarks need not have a very large 
splitting in this case. We have also shown that (in either of the methods used) the phase 
transition is not necessarily suppressed by large values of At and yU, which is encouraging 
since they are likely to be the sources of CP violation if electroweak baryogenesis indeed 
occurs in the MSSM. 

It seems clear that the two-loop corrections are more relevant for accurate results in the 
MSSM than for the standard model. We leave it for future investigation to see how the 
differences between DR and EP are affected when one includes the most important two-loop 
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corrections in DR, some of which have recently been incorporated into EP in references |]T2 



and |T^. A further shortcoming in the present state of calculations is that the DR method 
is no longer applicable in the MSSM if the fi or At parameters should be larger than several 
hundred GeV, or if the top squark becomes too light during the phase transition. While 
the latter could eventually be accounted for in 3D lattice simulations by including the light 
squark field in the effective action, to explore the possibihty of very large At or fi would 
require a simulation with the full 4D-theory. Let us finally point out, in favour of the DR 
approach, that for reliable computation of the phase transition dynamics one needs to know 
not just the order parameter v{T)/T, but also other quanitites like the latent heat and 
surface tension, which are not well-approximated by the EP approach in the standard model 
i- 
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